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Abstract 

The standard bounce formahsm for calculating the decay rate of a metastable 
vacuum cannot be applied to theories in which the symmetry breaking is due to 
radiative corrections, because in such theories the tree-level action has no bounce 
solutions. In this paper I derive a modified formalism to deal with such cases. As 
in the usual case, the bubble nucleation rate may be written in the form Ae~^. 
To leading approximation, B is the bounce action obtained by replacing the tree- 
level potential by the leading one-loop approximation to the effective potential, in 
agreement with the generally adopted ad hoc remedy. The next correction to B 
(which is proportional to an inverse power of a small coupling) is given in terms of 
the next-to-leading term in the effective potential and the leading correction to the 
two- derivative term in the effective action. The corrections beyond these (which 
may be included in the prefactor) do not have simple expressions in terms of the 
effective potential and the other functions in the effective action. In particular, the 
scalar-loop terms which give an imaginary part to the effective potential do not 
explicitly appear; the corresponding effects are included in a functional determinant 
which gives a manifestly real result for the nucleation rate. 
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1. Introduction 



When studying quantum field theories in a cosmological context one often en- 
counters situations where the field is not at the absolute minimum of the potential 
(the "true vacuum"), but is instead at some local minimum (a "false vacuum") 
which is higher in energy. The transition to the true vacuum state proceeds by a 
quantum mechanical tunneling process in which the field in a small region of space 
tunnels through the potential energy barrier to form a bubble of true vacuum. 
Once nucleated, the bubble expands, converting false vacuum to true as it does so. 

The bubble nucleation rate per unit volume, F, can be calculated by a method, 
due to Coleman,'^' which is based on finding a "bounce" solution to the classical 
Euclidean field equations. Thus, for a theory with a single scalar field, one must 
solve 



subject to the boundary condition that approach its false vacuum value as any 
of the x„ tend to ±00. F can then be written in the form 



where B is the Euclidean action of the bounce, while A is given by an expression 
involving functional determinants. In practice, the latter generally turns out to be 
quite difficult to evaluate, although one can show that A is equal to a numerical 
factor of order unity times a dimensionful quantity determined by the characteristic 
mass scales of the theory. 

A difficulty arises when one deals with theories where the symmetry breaking is 
a result of radiative corrections.'^' In such cases the true vacuum is not determined 
by V{(f)), but instead can be found only by examining the effective potential, V^fF(0). 
The bounce equation (1.1) is clearly inappropriate — in fact, if V^(0) has only a 
single minimum there will not be any bounce solution. An obvious alternative. 



(1.1) 



V = Ae 



B 
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which has been taken by a number of authors, is to modify this equation by 
substituting V^ff(0) for V{4>). 

Although plausible, and clearly a step in the right direction, this procedure 
raises some questions. The one-loop radiative corrections generate an effective ac- 
tion which contains not only V^g, but also terms involving derivatives (of all orders) 
of the fields. Can these terms be neglected when dealing with configurations, such 
as the bounce solution, which are not constant in (Euclidean) space-time? Even 
if this can be done in a first approximation, what are the nature and magnitude 
of the corrections which these terms generate? There are also questions relating 
to V^eff itself. First, the effective potential obtained by perturbative calculations 
differs considerably from that defined by a Legendre transform (the latter must be 
convex, while the former is not). The latter clearly does not lead to an appropriate 
bounce, but how precisely does the formalism pick out the former? Further, the 
perturbative effective potential is known to be complex for certain values of the 
fields. How is the imaginary part of V^fj to be handled? 

In order to answer these questions I develop in this paper a scheme for calcu- 
lating the bubble nucleation rate in theories with symmetry breaking by radiative 
corrections. The general idea is to use the path integral approach of Callan and 
Colemanj*' but to integrate out certain fields at the outset. This leads to a modified 
effective action which, although it differs somewhat from the usual perturbative 
SqS, gives a correct description of the vacuum structure of the theory and has a 
bounce solution which can provide the basis for a tunneling calculation. The fi- 
nal result for F may again be written in the form of Eq. (1.2), with A either of 
order unity or proportional to a small inverse power of the coupling. As might 
be expected, the leading approximation to B is obtained simply by replacing V 
with the dominant one-loop part of V^g in the standard procedure. The next-to- 
leading terms in V^f^, as well as the first corrections to the derivative terms in the 
effective action, then come in and give calculable and significant (i.e., larger than 
order unity) corrections to B. However, the corrections beyond this point do not 
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have a simple expression in terms of V^fj and the other functions entering S'cfj . In 
particular, the potentially complex parts of the effective potential do not enter the 
calculation directly, but only as part of more complicated functional determinants 
which can easily be shown to be real. 

This approach can be understood by recalling the Born-Oppenheimer approx- 
imation for the calculation of molecular energy levels. Because the natural time 

scale for the electrons is much shorter than that for the nuclei, the electrons can 
be treated as adapting almost instantaneously to changes in the positions of the 
nuclei. This leads to an effective action, involving only the nuclear coordinates, 
which may be used to describe situations in which the nuclei do in fact move slowly; 
thus, it can be applied to the vibrational spectrum of the molecule, but should not 
be used when studying the scattering of high energy particles off the nucleus. 

A similar situation arises in theories where radiative corrections change the 
vacuum structure. For definiteness, consider the case of scalar electrodynamics, 
where radiative corrections from photon-loop diagrams lead to spontaneous sym- 
metry breaking if the scalar field self-coupling is O(e^) and the scalar mass term is 
sufficiently small. This relationship between the couplings implies that the natural 
time scale for the variation of the electromagnetic field is much less than that for 
the scalar field. Consequently, if we are only interested in the long wavelength 
modes of the latter, as is the case when studying vacuum tunneling, we can in- 
tegrate out the electromagnetic field to obtain an effective action for the scalar 
field. 

The remainder of this paper is organized as follows. Sec. II reviews the method 
of Callan and Coleman. In Sec. Ill, I show how this method can be adapted to the- 
ories with radiative symmetry breaking. The treatment here is somewhat formal, 
so that it can be apphed to a variety of theories. The concrete implementation of 
the method in specific models is developed in Sees. IV and V, which discuss a sim- 
ple model with two scalar fields and scalar electrodynamics, respectively. Sec. VI 
contains some concluding remarks. 
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2. The Callan-Coleman Formalism 



Callan and Coleman'"' evaluated F by calculating the imaginary part of the 
energy of the false vacuum. This can be obtained from the quantity 

G{T) = (0(x) = (/)f,| e-^^ |0(x) = 0fv) 



(2.1) 



where the Euclidean action 



r/2 

S'(0) = J d^x J dx4 



(2.2) 



-T/2 



is expressed in terms of renormalized fields and parameters, while contains 
the counterterms needed to make the theory finite. Although divergent, these 
counterterms are of higher order in the couphng constants and, as usual, are treated 
as perturbations. 

The path integral is over all configurations such that takes its false vacuum 
value 0fy at ±T/2 and at spatial infinity. In the limit T — > oo Eq. (2.1) is dominated 
by the lowest energy state with a non- vanishing contribution (i.e., the false vacuum) 
and is of the form 

G{T) ^ e"^™ (2.3) 

where fl is the volume of space and S may be interpreted as the energy density 
of the false vacuum state. Because this is an unstable state, S is complex, with 
its imaginary part giving the decay rate, which in this case is simply the bubble 
nucleation rate. Dividing by fl gives the nucleation rate per unit volume, 

r = -2Im£ (2.4) 

The path integral may be approximated as the sum of the contributions about 
all of the stationary (or quasi-stationary) points of the Euclidean action 5'(0): the 
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pure false vacuum, the bounce solution with all possible locations in Euclidean 
space-time, and all multibounce configurations. In each case the contribution to 
the path integral is obtained by expanding the field about the classical solution 
(j){x): 

(f){x) = ^{x) + r]{x) (2.5) 

and then integrating over 77. To leading approximation one keeps only the terms in 
the action which are quadratic in 77. Expanding these in terms of the normal modes 
of S"{4>) — + ^"{4>) gives a product of Gaussian integrals. The terms of cubic 
and higher order in 77 can then be treated as perturbations, and have the effect of 
multiplying the Gaussian approximation for F by a power series in the coupling 
constants, with the first term being unity. The counterterms are also treated as a 
perturbation, with the leading contribution being from Sct{4>). 

Thus, the integration about the trivial solution (f){x) — 0fv yields 



Go = NKo e-^(<^fv)g-5ct(0fv)(i + . . .) 



(2.6) 



where N is a normalization factor. 



i^o = [det5"(0f,)]-V2 

= [det(-n + y"((/)f,))]-V2 



(2.7) 



and the dots represent terms, due to higher order perturbative effects, that can 
be neglected at the order to which we are working. Although the determinant is 
actually divergent, the divergences are cancelled by the factor of e"'^'^*. 

The calculation of the contribution from the bounce solution (phi^) is similar, 
but must be modified to take into account the fact that ^"'(^b) has one negative 
eigenvalue and four zero eigenvalues; for none of the corresponding modes is the 
integral truly Gaussian. The zero-frequency modes are treated by introducing 
collective coordinates corresponding to the position of the bounce. Integrating 
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over these gives a factor of QT and introduces a Jacobean factor which, for 0(4)- 
symmetric bounces, is given by 



J 



4^2 



(2.8) 



The negative mode is handled by deforming the contour of integration. Aside 
from a factor of 1/2, this leads to a contribution whose imaginary part is just 
that which would have been obtained from a naive application of the Gaussian 
integration formulas; i.e., 

ImG'i = G'oGb (2.9) 



where 



with 



and 



Gb = KiJnTe-''e 



ctl<Pb)— Oct(<Pfv) 



;i + ---) 



(2.10) 



(2.11) 



det'[-n + y"(0b)] 



-1/2 



(2.12) 



2 det[-n + y"(0fv)] 

Here det' indicates that the translational zero-frequency modes are to be omitted 
when evaluating the determinant. As before, the divergences in the determinant 
factors are cancelled by the terms containing S'ct- 

Finally, the n-bounce quasi-stationary points give a contribution of the form 
Gn — GoG^/nl. The bounce contributions then exponentiate, and one finds that 

-1/2 (2.13) 



S2 



-B 



47r2 



dct'[-n + l^"(0b)] 



dct[-n + V"{(f)fy)] 



-[<S'ct(0b)-<S'ct(0fv)] Q 



+ 



For comparison with later results we need to know the order of magnitude of 
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the various terms in this expression. To be specific, let us assume that we can 
identify a small coupling A and a dimensionful quantity a such that the potential 
can be written in the form 

= A(7^C/(V') (2.14) 

where U involves no small couplings and the dimensionless field ip — 0/(7. The 
minima of the potential must then be located either at = or at values of of 
order a. (With a single scalar field the most general renormahzable potential is of 
the form 

T/(0) = + + (2.15) 

The above assumption is then equivalent to assuming that c is of the order of A o", 
with a = m/y/X. ) 

By defining a dimensionless variable s — VXax, we may write the field equa- 
tions as 

(2.16) 

Prom the assumptions made above, this equation involves no small parameters 
and so has a bounce solution in which is of order unity and differs from the false 
vacuum within a region of a spatial extent (measured in terms of s) which is also 
of order unity. In terms of the original variables, the bounce has of order a and 
extends over a range of x of order l/(\/Xa). 

With the same change of variables, the action becomes 

Since the integrand contains no small parameters, while the volume of the bounce 
restricts the integration to a region of order unity, the bounce action B is of order 
A-i. 
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Similarly, the determinant factor Ki becomes 



Ki = - 



det'[(Aa2)(-a + f/"(^b))] 



-1/2 



det[(Aa2)(-n, + [/"(^f,))] 



-A a 



det[-n, + t/"(^fv)] 



(2.18) 



where the explicit factor of A^o"^ on the second line arises because the det' fac- 
tor involves four fewer modes than the det factor. With this factor extracted, 
the ratio of determinants is formally of order unity, although divergent.^ Finally, 
the Jacobean factor J is proportional to ~ A"^. Putting all of these factors 
together, we see that the nucleation rate is of the form 



4 ^C2/A 



(2.19) 



with ci and C2 both of order unity. 



1 It was asserted above that these divergences are cancelled by divergences in Set (^) — Set i4>fv ) • 
To see that these counterterms are of the correct order of magnitude, note that with our 
assumptions the one-loop contributions to the counterterm Lagrangian are of order X^a^; 
since the bounce has a spatial volume of order l/(-\/Acr)'*, the difference in the counterterm 
actions is of order unity, which is what is required. 
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3. A Formalism for Theories with Radiative Symmetry Breaking 

I now describe how the methods of the previous section can be adapted to 
theories where radiative corrections to the tree-level potential qualitatively change 
the vacuum structure. The general formalism is derived in this section, while the 
detailed implementation for two specific theories is described in Sees. IV and V. 

In theories of this sort the fields may be divided into two sets, denoted Xj 
The former give rise to the radiative corrections which effect the changes in 
the vacuum structure, while the various vacua are distinguished by the value of the 
latter. In order that the radiative corrections from the x-loops be comparable to 
the tree-level terms, the interactions among the 0j must be weak compared to the 
(f)-X interactions. Hence, the Xj ^lay be viewed as "fast" degrees of freedom which 
can be integrated out in a Born-Oppenheimer type approximation. For simplicity 
of notation, I will assume here a theory with only field of each type, both assumed 
to be scalars; the generalization to more complicated theories is straightforward. 

As usual, the bubble nucleation rate is obtained from the quantity 

G{T)^ Jmdx]e-'^'l'^^^ (3.1) 

(To simplify the equations the counterterm action has not been explicitly written.) 
Because S{(f), x) does not display the correct vacuum structure, its classical field 
equations do not have a bounce solution, and so the methods of Sec. II cannot 
be directly applied. This difficulty can be circumvented by integrating over x to 
obtain 

G{T) = J [d^Je-^^"^) (3.2) 

where 

W{(f))^-]n j [dx]e-^^'t>,x) (3.3) 

may be thought of as an effective action. W has a simple graphical interpretation. 
The integral on the right hand side of Eq. (3.3) is equal to the sum of all vacuum 
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Feynman diagrams in a theory of a quantized x-field interacting with a fixed c- 
number 0-field; i.e., the sum of all graphs with only internal ^-lines and external 0- 
lines. Its logarithm gives the sum of all connected vacuum graphs. In general these 
include one-particle reducible graphs, although in many cases these are eliminated 
by symmetry considerations. 

Vr(0) should be compared with the more familiar effective action, S(;ff{(f),x), 
which generates one-particle irreducible Green's functions. For theories in which 
W{(f)) only receives contributions from one-particle irreducible graphs, Scs{(f),0) 
and W{(j)) differ by the contribution from graphs with internal 0-lines. Our as- 
sumptions about the size of the 4> self-couplings imply that these graphs are sup- 
pressed. Hence, the dominant terms in 5'eff(0, 0) and W{(f)) agree, and so the latter 
correctly reflects the vacuum structure of the theory. 

One might therefore envision using the classical field equations implied by W 
to determine a bounce solution which would be the basis for a nucleation rate 
calculation. Two practical difficulties arise. First, one cannot in general obtain 
a closed form expression for W, but only a perturbative series. Second, W{4>) is 
a nonlocal functional, so that even if a closed-form expression were available, the 
equations determining its stationary points would be rather unpleasant. 

One way around these difficulties is to find a local action VFo(0) which is a 
sufficiently close approximation to W and then use its stationary points as the 
basis for the nucleation rate calculation. The possibility of such an approximation 
arises from the fact that the nonlocality of W is significant only on distances 
shorter than a characteristic size set by the x-field interactions. For 0-fields which 
are slowly varying relative to these scales, W{(f)) can be expanded in a "derivative 
expansion" of the form 



where the dots represent terms with four or more derivatives. (A similar expansion 
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(3.4) 



of the effective action, 



(3.5) 



is often made. In line with the previous remarks, t^(0) and ^(0) differ from 
^ff (0) 0) and Z(j){(f), 0) by the omission of graphs with internal 0-lines.) 

The desired approximate action is obtained by keeping only the first two terms 
in the derivative expansion of W^(0), and then only the lowest order contributions 
to these; thus 

~1 



W^o(0) 



d X 



;(5^0)' + ^l-loop(0) 



(3.6) 



Here Vi 

-loop(0) is the sum of the tree-level 1^(0) and the contributions from graphs 
with a single x-loop; to this order it is equal to 14ff(0, 0). Since the loop correc- 
tions included in Vi_ioop(0) are precisely those responsible for altering the vacuum 
structure of the theory, Wq possesses an appropriate bounce solution. 



Thus, for slowly varying we may write 



W{(f))^Wo{(f))+SW{(l)) (3.7) 

with 5W{(I)) representing subdominant terms. We can (and will) make the same 
decomposition for arbitrary 0, but 5W will then not necessarily be small. In 
particular, because the path integral includes contributions from fields with Fourier 

components of arbitrarily high momentum, for which the derivative expansion is 
not valid, 5W cannot be treated as a small perturbation within the path integral. 

Nevertheless, the path integral may be approximated by expanding about the 
stationary points of Wq, provided that these are themselves slowly varying func- 
tions. Thus, let us suppose that the field equation implied by Wq, 

□0 = ?%J22P (3.8) 
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has a solution (f){x). Expanding W about this solution gives 



W{(l)) = W{(l)) + J d^z W\(t)] z)r]{z) + ]^jd^z d^z' W"{4)- z, z')r]{z)ri{z') + 0{ri^) 

(3.9) 



where 'r]{x) — (j){x) — (f>{x), and primes denote variational derivatives; e.g., 

SW{(j)) 



W'{^;z) 



(3.10) 



(The term linear in rj is nonvanishing because is not a stationary point of the 
full W{(l)).) 

We now insert Eq. (3.9) into the path integral. The cubic and higher order 
terms in 77 can be treated as small perturbations. If these are dropped, the integral 
becomes Gaussian. About the false vacuum, (j){x) = 0fv, W"{(j)) has no zero or 
negative eigenvalues and the integration gives 

Go = e-^^t>^^ ^w'{ct>^)[W"{<t>^)rw'{4>,.) ^^^^ iy"(0f^)]-i/2 (3.11) 

(For the sake of compactness, the spatial arguments of the variational derivatives 
and the associated spatial integrations have not been explicitly displayed.) As in 



the standard case, the bounce solution has zero frequency modes which must be 
handled 
leads to 



2 

handled separately. Introducing collective coordinates and proceeding as usual 



ImGi = lQTe-^(<^'') eW^'(^Or"(0b)]-T^'(<^b) |detV'(0b)r'/' J (3.12) 

where now det' indicates that the determinant is restricted to the subspace or- 
thogonal to the zero modes of W^. Similarly, [yV")^^ is to be evaluated in this 



2 Because the bounce is obtained from the truncated action Wo, the usual translational modes, 
proportional to 9^(/)b, differ slightly from the zero modes of W" \ however, it can be shown 
that the effects due to the difference are suppressed by a power of the coupling constant 
and hence can be neglected at the order to which we will be working. 
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subspace, while the Jacobean J is given by 



The path from these results to the nucleation rate is just as before, and gives 



detiy"(0fv) 



J(l + ---) (3.14) 



where 



Ci = W{4>^) - W{4>^) (3.15) 



and 



C2 = W\cj>^)[W"{cl>^)]-^W\cj>^) - {ci>^ ^ 0f,) (3.16) 



The next step is to expand the various terms in Eq. (3.14) in powers of the 
coupling constants. The details of this depend on the structure of the theory and 

the magnitudes of the various couplings. The general features described below hold 
for the specific theories considered in the next two sections and are likely to be 
true for all theories with radiative symmetry breaking. As in Sec. 11, the expansion 
of r is only carried out to terms of order unity; terms which do not contribute to 
that order are represented by dots. 

First, the bounce solution has a characteristic length scale which is large com- 
pared to those characterizing the X'X ^iid interactions. It therefore varies 
slowly enough to allow a derivative expansion of W{(f)]^). (We will see that in some 
situations this expansion must be terminated after the first few terms, leaving a 
nonlocal remainder. This does not materially affect the picture outlined here.) For 
both = and = 0fv, it is then possible to write VF(0) as a sum of terms of 
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successively higher order in the couphngs: 

W{^) = Wo{^) + Wi{^) + W2{^) + ■■■ (3.17) 

with Wq given by Eq. (3.6) and the remaining terms containing higher order contri- 
butions to the derivative expansion. For example, in the theories considered below 
Wi contains two-loop contributions to V and one-loop contributions to Z. 

A similar expansion can be used for W'{(f); z). The zeroth order term vanishes 

because is a stationary point of Wq, and we have 

W\4>\z) = W[{4>\z) + --- (3.18) 

Matters are less simple for W"{4); z, z'). Although the derivative expansion can 
be used when \z — z'\ is large, the behavior for small \z — z'\ is sensitive to the 
high momentum modes and so the derivative expansion fails even for constant 0. 
However, the relation 

W" = [1 + {WHy^SW"] (3.19) 

can be used to obtain formal expansions for {W")~^ and detW" as power series 
in (Wq)^^5W" . The actual utility of these expansions depends on the size of 
the contribution from the region of small \z — z'\. In the calculation of C2 this 
contribution is subdominant and 

C2 = W^l(0b) W(0b)]"'W^l(0b) - (0b ^ 0fv) + • • • (3.20) 

For the determinant factor, on the other hand, more terms must be retained: 

det'[W^"] = det'[V^,^'] det[/ + {Wi;)-^6W"] 

= det'[W^'] exp {tr ln[7 + {WHy'^SW"] } ^3 2i) 

= det'K'] exp Lt{WI;)-^6W" + ^tr [{WHy^dW'f + ■ ■ - j 
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Inserting these expansions into Eq. (3.14) leads to 



dctV(^^(0b) 
detiy^'(0fv) 



1/2 




r = e~^° e~^' e~ 



47r2 



(! + •••) 



(3.22) 



where 



Bo = W^o(0b) - W^o(0fv) 



(3.23) 



Si = m(0b) + ^tr[<(0b)] 



'5W"{<I>^) - (0b ^ 0fv) 



(3.24) 



and 



^2 = W-2(0b) - W-o(</'b) K'(</'b)]-'W-o(</'b) + ^tr { W(</-b)]-'(5iy"(0b)}' 



- (0b 0fv) 



(3.25) 



We thus have the naive expression for F, obtained by the standard procedure 
with y(0) replaced by the one-loop effective potential, multiplied by a correction 



couplings, while B2 is of order unity. Thus, although the B2 correction is of the 

same order of magnitude as the determinant factor, which in practice cannot be 
evaluated precisely, the Bi factor modifies the naive result in a significant and (at 
least numerically) calculable manner. 

The form of these correction factors can be understood. Bq is obtained by 
evaluating the approximate action Wq at its stationary points, whereas what we 
really want is the exact action W evaluated at its own stationary points. The 

difference between these two quantities can be obtained by treating 6W = Wi + 
W2 + ■ ■ ■ as a small perturbation; up to second order, this gives the first term in Bi 
and the first two terms in i?2- As we will see below, the remaining terms correct for 
the fact that W is not the full effective action and take into account contributions 
to the full effective action which would otherwise be neglected. 



factor e (^1+^2). We will find that Bi is proportional to an inverse power of the 
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4. Example 1: A Scalar Field Theory 



As a first example, consider a theory with two scalar fields governed by the 
Lagrangian 



with 



(4.2) 



Both jj?' and are positive. In order that the x-loop corrections to the efi^ective 
potential be able to generate a symmetry-breaking minimum at a value (f) — a 0, 
A must be while / should be 0{g^) to ensure that the x ~^ ~X symmetry 

remains unbroken. There is also an upper bound of order g^a on which will be 
displayed in detail below. Finally, it will be convenient to assume that m is 0{ga). 

Renormalization conditions must be specified in order to fix the counterterms 
in Cct- For the lowest order approximation to the effective potential only the 
definitions of fj? and A are needed. These can be given by choosing an arbitrary 
renormalization scale a and requiring that 



d^Vr 



cff 



(4.3) 



=x=o 



and 



A = 



d^Vr. 



eff 



(4.4) 



For X = the effective potential is given by the sum of graphs with only exter- 
nal </)-lines, all carrying zero four-momentum. The leading contributions are from 
tree graphs and graphs with one x-loop, as well as the corresponding counterterms. 
Let us denote the sum of these by V^i-ioop(0)) although it should be kept in mind 



3 Since x = is to be understood for the remainder of this discussion, T4ff and related 
quantities will be written as functions of a single variable. 
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that this does not include the graphs with one 0-loop. We then have, up to an 
additive constant, 



In 



1 



^ ) + 



1 /X N 1 //^ 



0^ 



(4.5) 



where 



and 



1/1 



In 



^1 



Stt^ 



167r2 



2-2 



12g'a 



(4.6) 



(4.7) 



It is most convenient to choose the renormahzation point a to be u, the location 
of the minimum of Vi_ioop- The requirement ^i-ioopl*^) ~ ^ then relates A to the 
other parameters via 



A = -3 



+ yi 



and yields 



(4.8) 



14-loop(0) 



647r^ 



/ 2 I 2 /2\2 

(m + ) 



In 



+ g'^cj)^ 
+ g^a^ 



1 //^ ^ 



(4.9) 



For (T to be a minimum, ^i^-ioopl*^) i^^st be positive, which leads to the bound 



2 , 9'<y' 



4- 



m 

gr2(j2 



In 



(4.10) 



A slightly stronger bound on // is obtained if one requires, as we will, that the 
minimum at = cr be deeper than that at = 0, so that the latter corresponds to 
a metastable false vacuum. 



17 



The higher order contributions to the effective potential can be ol:)tained from 
Feynman graphs with "dressed'' propagators that take into account the bacl<;ground 
0-field. Thus, the effect of summing over all numbers of interactions with a constant 
external 0-field is to replace the standard propagator (/c^ — m^)~^ by 

[k^ -rr?- g^ct)'^]-^ = [k'^ - M^{(f))]-^ (4.11) 

Because the x-loop contributions to \4ff are the of the same magnitude as the tree 
terms, the dressed ^-propagator is somewhat more complicated. In order that the 
perturbative order of a graph increase with the number of loops, these one-loop 
terms must be included in the dressed propagator, which then takes the form 

[k' - n-ioop('/')]"' ^ [k' - M\<pr' (4.12) 

This inclusion of loops in the propagator means that the loop-counting for a graph 
with dressed propagators will not necessarily agree with that for the corresponding 
graphs with elementary propagators. Thus, the one 0-loop contribution to the 
effective potential will include graphs which are multiloop when drawn with ele- 
mentary propagators. Some examples of these are shown in Fig. 1. In the last of 
these the momentum assignments have been shown explicitly in order to emphasize 
that they are not the ones usually associated with this graph. It should be clear 
from this example that graphs with both x ^-nd internal lines must be analyzed 
carefully to avoid double-counting. 

The next contributions to V^fj, of order fg'^ and g^, arise from the two-loop 
graphs in Fig. 2, the one-loop counterterm graphs in Fig. 3, and the 0{fg'^) and 
0{g^) counterterms. The double-counting issue arises with the second graph of 
Fig. 2, which has the same topology as the graph of Fig. Ic, but a different (i.e., 
the standard) assigment of momenta. However, since the latter graph is 0{g^) 
(after renormalization) , the issue can be ignored at this point. 
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The sum of these graphs gives a contribution which, in terms of an arbitrary 
renormahzation scale a, can be written in the form 



204^? [/(^' + 9V) + 32/02] + ^2^2) 1^2 



+ 



2^2m^2 (rif±9^ 



{cif + C2g'^){m'^ + g 



647r2 

+ (C3/ + 04^)^2 (m2+/02 



In 



2j,2 



In 



(4.13) 



777.2 _|_ g2g.2 



4 



Here the Ck are numbers of order unity which depend on the precise choice of the 
renormahzation conditions fixing m?, g^, and the x-field wave function renormal- 
ization, while y2 and Z2 are to be chosen so as to satisfy Eqs. (4.3) and (4.4). If a 
is again chosen to be a, with the latter now understood to be the minimum of the 
two-loop approximation to V^fj, this result combines with li_ioop to give 



^2-ioop(0) = ^(1 + cif + C2g')(m' + g^'^f 



In 



m 



2a2 



m? + g'^cj'^ 



2j.2 



+ (cs/ + C4,g )m {m +g 



2i2\ 



In 



^m2 + 5'2(j2 



m + g 
w? + g'2cr2 



(4.14) 



where y^yi + y2- 

Among the contributions at 0{g^) are those from the graphs with a single 
0-loop. These include a term proportional to ^'^yi^\oo^{'P) which is complex when- 
ever V/^ 

-loop negative. While such complex effective potentials are not unfamiliar, 
and their physical interpretation is understood,'^' their implications for bubble nu- 
cleation calculations have remained somewhat unclear. We will see below how the 
potential problems associated with these are avoided. 
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Let us now turn to the bubble nucleation problem. The first step is to calculate 
W{(j)). Because the Lagrangian only contains even powers of x, only one-particle 
irreducible graphs contribute. One obtains 



+ counterterms + • • • 

(4.15) 

where the contributions from tree and one-loop graphs have been shown explicitly. 
The second term, arising from the loop graphs, is sensitive to the value of 0(y) 
over a region of size ~ M~^{(j)) ~ {ga)~^. Hence, for our purposes a slowly varying 
(j) will be one which varies slowly over this distance. Expanding such a field about 
its value at a point x leads to the derivative expansion for W. 

The first term in the expansion, V, is obtained by treating as a constant. To 
order g"^ this is the same as the standard calculation of the effective potential, and 
gives the result displayed in Eq. (4.9). However, the two quantities differ at order 
g^, since the two-loop graph of Fig. 2b contributes to V^g but not to V. Because 
this graph is divergent, while the counterterm contributions to the two potentials 
are the same, V, unlike V^g, cannot be finite. This divergence is acceptable because 
V does not directly correspond to any physical quantity. Indeed, it combines with 
other divergent quantities to give a finite result for F. 

The calculation of the next term in the derivative expansion is equivalent to 

extracting the term of order from the sum of graphs with two external 0-lines 
carrying nonzero momentum p and any number of zero-momentum external 0- 
lines. Using dressed propagators, the leading part of this sum can be represented 
by the single graph of Fig. 4, whose evaluation gives 



fip') = /(O) + 



87r2 



dx In 



1 + 



p'^x(l — a 



(4.16) 
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For ^ M'^^cJ)) this can be expanded in powers of p^: 



/(pV/(0) + j|^ 



9 4 
P p 



M2(0) 5M4((/)) 



(4.17) 



(The failure of this expansion when <^ M'^{(f)) reflects the failure of the derivative 
expansion when the background field is not slowly varying.) The first term, /(O), 
has been included in the calculation of V. The 0{p^) term leads to 



ZU) = 1 + ^ 



+ k 



+ 0{g^) (4.18) 



where is a constant of order unity which is chosen to enforce the wave func- 
tion renormalization condition. (Note that to this order ^(0) = ^,^(0,0).) The 
O(p^) term gives a contribution to the derivative expansion of the form H {(/)) (n(p)'^ ■ 
In addition, there are other four-derivative terms, proportional to ((9^0)^ 00 and 
(9^0)2 (5,y0)2, which can obtained from graphs with three or four external lines 
carrying nonzero momentum. 

The bounce solution is determined by the approximate action of Eq. (3.6) 
with Vi-ioop(0) being just the V^i_ioop(0) of EQ- (4-9). Scaling arguments similar 
to those of Sec. 11 show that it has a spatial extent of order {g'^a)~^ and thus 
varies slowly enough to allow a derivative expansion. By contrast, in a theory 
with couplings of "normal" magnitude (A ~ and /j? ~ fi'^c"^) the bounce would 
have a size ~ {go')~^ and would not be slowly varying. Of course, in that case the 
loop corrections to the effective potential would not change the vacuum structure, 
so the standard calculation of F could be used and there would be no need for a 
derivative expansion. 

We need to know the magnitude of the various terms in the derivative expan- 
sion of iy(0b). Since 0b ~ cr, we see that V{4>h) ~ 9^<^^j while our estimate of the 
bounce size implies that (9^0b)^ ~ More generally, consider a term contain- 

ing j derivatives and k explicit factors of 0b, which might be written schematically 
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as 

= HjMd^4>t (4.19) 

Hjif((j)) must have dimensions of {mass)^~^~^; since it is derived from graphs with 
only internal lines, the mass appearing here must be the effective x mass, M{4>). 
Further, the structure of the graphs imphes that there is a factor of g for each 
explicit factor of (p, except in the tree-level contribution to Z. Finally, each deriva- 
tive brings in a factor of g^a. Thus, apart from the terms contributing to Z (i.e., 
the case j — k — 2), 

where the fact that M'^{(j)) > ~ g^a"^ has been used and it has been assumed 
that (f) ^ a. 

We can now estimate the various terms entering the formula for F, beginning 
with the quantities Wj in Eq. (3.17). Because F involves the differences between 
these quantities evaluated at = 0b ^-nd at = 0fv, the integrals are effectively 
restricted to a region of volume ~ (l/gf^cr)^. Using the above estimates for the 
integrands, we then find that the first two terms are 

Wo^ J d^x[Vg4 + (90)2] ~ g-^ (4.21) 

and 

Wi = J d^x[Vg6 + Zg2{d4)f] - (4.22) 

Here the subscripts indicate the order of the contributions to V and Z, with Vgi 
being the same as Vi-ioopj ^-nd V^e including both and fg^ terms. The next 
term, IF2, is of order unity and contains Vgs (arising from graphs with up to three 
loops), iZgi (arising from one- and two-loop graphs) and the one-loop contributions 
to the four-derivative terms. 
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We also need expansions of the first two variational derivatives of W. In 
principle, the variations should be performed before the expansion is carried out, 
although this turns out not to matter for the first variation. Thus, to one-loop 
order 

W'{(l); z) = -\D(t){z) + V'{(t){z)) + g^(t){z){z\ [-□ + m2(0)] ~^ \z) + counterterms 

(4.23) 

The leading contribution to the third, nonlocal, term is obtained by taking to be 
constant. This combines with V' to give V[_Yoop^ ^^^^ leading approximation 
to W' is indeed the variation of the leading approximation to W . One can check 
that analogous statements hold at the next order, thus verifying Eq. (3.18). 

Up to one-loop order the second variational derivative is 

W"{<t>; z, z') = 5(4) (^ - z')[-U + V"m + ^2^(4) (^ - z'){z\ [-□ + m2(0)] \z) 
-2g^(p{z)(p{z ){z\ [-n + M'^{(j))Y^ I z') 2 counterterms 

(4.24) 

For \z — z'\ ^ M(0) the nonlocality of the third term cannot be ignored and the 
derivative expansion fails, even in a spatially constant background. To see this, let 
us suppose that (t){x) is constant, so that we can work in momentum space with p 
being the momentum conjugate to z — z'. The last term in Eq. (4.24) is then given 
by the graph of Fig. 4. If p were set equal to in the evaluation of the loop, then 
the one-loop terms shown here would give Wq, which in momentum space is simply 
Wq{p) ^p^ -\- V^^ioop(^) = + M^{(p). Subtracting this and recalling Eq. (4.16), 
we see that 5W"{p) ~ g^Vlnfl -Fp2/M2(0)]. If p^/M'^{(t)) is small, this can be 
treated small perturbation, but this is not the case when p^/M^(0) is so 

large that the logarithm overcomes the factor of g^. Correspondingly, in position 
space (2;|[l^Q']~^5iy"|2;') is small except in a region \z — z'\ <^ M~^{(f)). Similar 
estimates clearly apply with a slowly varying background field. 

W" enters the formula for the nucleation rate both through the quantity C2 
defined by Eq. (3.16) and through the determinant factor. In both cases Eq. (3.19) 
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leads to a formal expansion in which the higher order terms are suppressed if 
6W" is indeed a small perturbation; i.e., if the contribution of the short-distance, 
high-momentum region is suppressed. Roughly speaking, C2 can be viewed as 
corresponding to a tree graph in which two factors of SW' are connected by a 
propagator. The magnitude of the momentum running through this propagator is 
set by the spatial extent of the bounce solution, and is indeed small enough for SW" 
to be treated as a perturbation. (This can be checked by a detailed examination 
of the contribution from the small z — z' region.) We may therefore write 



C2= / d4zdViy^(0b;^)w^l(0b;^')(^lW(</'b)]"'|^')-(</'b^0fv) + --- (4.25) 



The size of the bounce restricts the z integration to a volume of size 
while the falloff of [H^q]^ restricts z — z' to a, similar volume. The two factors 
of W{ are each of order g^a^, while [Wo']"-*^ gives a factor of roughly {z — z')'^. 
Combining these facts, one finds that the terms shown exphcitly in Eq. (4.25) are 
of order unity. 

For the determinant factor, we recall from Eq. (3.21) that det W" can be writ- 
ten as a product of det Wq and the exponential of a power series in 
The former combines with the Jacobean to give times a factor of order unity, 
as in the standard calculation. The latter corresponds to a sum of one-loop graphs 
in which a number of insertions of SW" are connected by an equal number of 0- 
propagators. This sum has an ultraviolet divergence which can be attributed to 
the graphs with one and two insertions of SW"; even after the cancellation of these 
divergences by counterterm contributions, we should expect the high momentum 
region to be enhanced, and these terms to be anomalously large. This is in fact 
what happens. Working in momentum space, and ignoring the spatial variation of 
(j) for the sake of clarity, we may write the nth term in the expansion as 



The momentum integration is quartically divergent; after cancellation of infinities 
by counterterms a finite term proportional to the fourth power of M(0), the largest 



tr 
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relevant mass, remains. There is a factor of for each 5W" . Finally, because we 
must eventually take the difference between the contributions of the bounce and 
of the pure false vacuum, the x integration is effectively restricted to a volume 
~ (gf^cr)"^. The net result is a contribution of order g'^^~^, implying that both the 
n — 1 and the n — 2 terms are of order unity or larger, as indicated in Eq. (3.21). 

The dominant, 0{g~'^), contribution from the n = 1 term is obtained by treat- 
ing as a constant. This term then essentially reproduces the two-loop graph of 
Fig. 2b, which contributes to Vefj but not to V. (The definition of Wq actually 
requires a subtraction corresponding to the graph of Fig. Ic, but this is simply the 
correction for double-counting noted previously.) To this order, then, the net effect 
is to simply replace V by the full effective potential. Furthermore, since Z and 
are equal at O(e^), it might seem that the derivative expansion of W{(f)) is simply 
being replaced by the derivative expansion of Ses{(f), x = 0). 

However, when we go to the next higher order it becomes clear that this is not 
the case. At 0{g^), Ves differs from V by the contribution from two- and three- 
loop graphs with both x ^-nd (p internal lines plus the contribution from graphs 
with a single 0-loop. The former set of graphs is reproduced by the n — 2 term in 
the expansion of the determinant and by the next to leading contribution of the 
n = 1 term. The contribution of the latter set is contained in the factors of det Wq'; 
indeed, for a constant background (p the logarithm of this determinant is just the 
spatial integral of this part of the effective potential. Matters are more compli- 
cated when this determinant is evaluated in the background of the bounce. Since 
the effective 4> mass is a factor of g smaller than the effective x mass, the deriva- 
tive expansion of det Wq requires a more slowly varying background field than is 
needed for the expansion of the terms in VF(0) arising from x-loops. Precisely 
because the bounce is determined by the field equations of Wq, its characteristic 
size is incompatible with a derivative expansion of det H^q (0b)) ^^'i so the 0-loop 
contribution to V^g cannot be simply isolated. 

This resolves a puzzle associated with the complexity of the perturbative ef- 
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fective potential. As noted above, the 0-loop contribution to K;ff(0) is complex for 
values of such that Vi_ioop{4>) is negative. Since the bounce solution 0b (^) li^s 
in this range for some values of x, a manifestly complex result for the nucleation 
rate would have been obtained if the calculation had explicitly involved the full 

4 

effective potential. 

Of course, the fact that this effective potential contribution cannot be isolated 
does not rule out the possibility that a complexity problem could arise with the full 
determinant, where it would be manifested by the existence of a negative eigenvalue 
of VF(3'(0b) in addition to the one which occurs for any bounce. However, from the 
fact that the bounce corresponds to the solution of a minimization problem (that 
of finding the path through configuration space with the smallest WKB tunneling 
exponent), it is easy to show'*^' that the bounce solution of minimum action can 
never have additional negative modes. Hence, jdet'l^Q (^b)^^^ is real and the 
imaginary part of the effective potential has no effect on the bubble nucleation 
calculation. 

To summarize these results, we may write 

r = ie-(-^°+^i) (4.27) 



where 



Bo 



- (0b ^ 0fv) ^ = 0{g 



(4.28) 



d^x 



Kff('^b) + ^^|'(0b)(5M'/'b)' 



(0b ^ 0fv) = 0{g-^) (4.29) 



and the pre-exponential factor A, now understood to include the contribution from 
the 0(1) part of the exponent in Eq. (3.22), is equal to cr^ times a dimensionless 
factor of order unity. 



4 Although det W,^'((/)fv) can be expanded in a derivative expansion, no complexity problem 
arises because Vea{(j)u) is real. 
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5. Example 2: Scalar Electrodynamics 

[^] 

A second example is provided by scalar electrodynamics. In terms of the real 
and imaginary parts, 4>i and 02 of the scalar field, the Lagrangian takes the form 



where 0^ = 02 _|_ ^2 order that the one-loop radiative corrections give rise 
to the desired vacuum structure with both a symmetric minimum at = and 
a symmetry- breaking minimum at = cr, we require that A be O(e^) and that 
iu2 be positive and 0{e^a'^). As with the previous example, fi and A are to be 
fixed by Eqs. (4.3) and (4.4), while the other renormalization conditions will not 
be explicitly given. 

It is most convenient to calculate the effective potential in Landau gauge. In 

this gauge all graphs with a zero-momentum external 0-line entering a 0-0- 

vertex vanish; this leads to a considerable reduction in the number of graphs to be 
5 

calculated. A further simplification follows from the observation that V^fj can only 
depend on 02, so that we may calculate it with 02 = and 0i = 0. To 0(6^) the 
result is 



4 
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(5.2) 



If, as before, we choose a — a, the minimum of Vi_ioop, we can eliminate A and, 
after adding a constant, write 



14-loop(0) = ^0 



1" 


I/.2 


2_ 


4(72 



(02 -(.2)2 (5.3) 

Requiring that l^i^ioop('^) > implies that ji'^ < Se^a'^ /Idn'^; for this asymmetric 



5 The issue of gauge dependence is an important one, with some issues still to be settled. 
It is known that the effective potential is gauge-dependent beyond lowest order; this is ac- 

[7] 

ceptable becauwc it is not in general a directly measurable quantity. Physically measurable 
quantities, on the other hand, should not depend on the choice of gauge. However, an inves- 
tigation of the gauge-dependence of the bubble nucleation rate in scalar electrodynamics 
found an apparent dependence on the gauge choice. 
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minimum to be the true vacuum, the stronger condition /j? < 3e'*(T^/327r^ must 
hold. 

As in the previous example, the higher order contributions to the effective 
potential are most easily obtained from Feyman graphs with modified propagators. 
In Landau gauge with 02 = these are the dressed photon propagator with an 
effective photon mass 



(5.4) 



and the dressed scalar propagators with masses given by 



loop 



dcl>i 



bi=4>,4>2=Q 



(5.5) 



and 



loop 



90 



i>i=4>,4>2=Q 



(5.6) 



In gauges other than Landau gauge there would also be mixed 02-^^ propagators. 
As before, the inclusion of loops in the dressed scalar propagators implies that 
caution must be used to avoid double-counting. 

The O(e^) contributions to V^eff come from the two-loop graphs of Fig. 5, one- 
loop counterterm graphs identical to those of Fig. 3, and the O(e^) counterterms. 
Combining these with the O(e^) terms, with the renormalization point now chosen 

|"9| 

to be the minimum of the two-loop approximation to the effective potential, gives 



^2-loop(0) 



3e^ 
647r2 



(l + ce2)0^ 



In 



a 



+ 



5e^ 



5127r' 



-0Mn2 



4(72 



(02 - a')' 



(5.7) 



where c, a number of order unity, depends on the precise specification of the renor- 
malization conditions. 
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The contributions from purely scalar loops enter at 0(e). These make the 
effective potential complex, since Aif is negative in the region where ^i^_ioop < 0, 
while M2 < when \4>\ < a. 

As in the previous example, the function Z(j,{(j)) will also enter the calculation 
of r. At O(e^) this receives contributions from the graphs of Fig. 6, and is given 

by 

3e2 

where of order unity, enforces the wave function renormalization condition. 

Let us now turn to the calculation of the bubble nucleation rate. Before going 
into the details, we can use the example of the previous section to anticipate a 
number of the results. The bounce is presumably determined by field equations 
involving Vi_ioop- It is then easy to see that it must have a constant phase, and 
so by a global gauge rotation can be made entirely real. Furthermore, the spatial 
extent of the bounce must be of order (e^cr)"^, thus making it slowly- varying 
enough to justify a derivative expansion of the photon-loop terms, but not of those 
due to graphs with purely scalar loops. This will make it possible for the O(e^) 
part of the effective potential to enter the calculation of F in a straightforward 
fashion, while still avoiding the troublesome 0{e^) terms. 

There are at least two possible approaches to the calculation. Since the sym- 
metry breaking arises from radiative corrections involving photon loops, the most 
straighforward procedure would be to integrate over the gauge field to obtain an 
action VF(0i, 02)- Since enters the Lagrangian at most quadratically, this inte- 
gration can be done exactly, yielding 



In 



0^ 



(5.8) 
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+ ^ tr y d^x {x\ In [- + + 1^) 

d^xd''x'f{x)f{x'){x\ [-5^,u + d^d^. + 5^,M\cl>)Y^\x') 
+ counterterms 

(5.9) 

where = e(02f^^0l — 0l<9^02) and the factors of S^i, arise from the Euchdean 
metric. The second term on the right hand side corresponds to graphs with a single 
photon loop, and is analogous to the one x-loop term in the model of Sec. IV. The 
last term corresponds to a graph with a photon line connecting two currents. While 
it vanishes for both the constant false vacuum solution and the bounce, it does not 
vanish in general and, indeed, it is clearly needed to bring the O(e^) part of the 
effective potential into the exponent of F. However, because it cannot be expanded 
as a sum of local terms, the analysis of this term introduces complications not 
encountered in the previous example. 

These complications are avoided by an alternative approach. Since the bounce 
can be chosen to be entirely real, it should be possible to integrate out 02, as well 
as Ajj,, from the very beginning to obtain an effective action 



d X 



^(a,0)^ + i;.V + T'^^ 



A 



+ -tr In G + counterterms H (5.10) 

2 



Here the dots represent multiloop contributions while G is the matrix 



G 



6^26 



SA^,S<I>2 



fj,2 



--A^=0 



\G2u{(t>) G22{4>) J 



(5.11) 



In Landau gauge with a constant background (f)i{x) the off-diagonal entries of this 
matrix vanish. For a nonconstant but slowly- varying background one can write 



G 



G^M) 

G22{(P) 



1 + 







G-}{cl>)G^2{cl>) 



G^i{<P)G2.{<P) 



(5.12) 
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and then expand W as 



I (fx 



+ ^ tr In G,j,i, + In det G22 + ^ tr G~^G 1^2022 (^2^ + counterterms H 

(5.13) 

The first (tree-level) term and the second term, corresponding to one-photon-loop 
effects, dominate. Aside from a minor point noted below, these can be treated just 
like the analogous terms for the previous model. Let us write these as Wq + SiW, 
with Wq, given by Eq. (3.6), again being the approximate action which determines 
the bounce solution (t)h{x). The spatial extent of the bounce is such that 5iH^(0b) 
can be expanded in a derivative expansion. This expansion does not give any 
further contribution to V, but does give the contribution to Z corresponding to 
the graph of Fig. 6a. 

Let us consider next the fourth term, which may be denoted 52W. This cor- 
responds to one-loop graphs with both a photon propagator and a 02-pi'opagator. 
Because G22 is obtained from the fundamental action S, the latter propagator will 
have a mass A^2(0) given by the tree- level potential rather than by Vi_ioop, as is 
the mass Ai2{4') of Eq. (5.6). The replacement Ai2 M.2 is obtained by sum- 
ming over multiloop graphs with arbitrary numbers of photon loops attached to 
the 02-line, as in Fig. 7; let us denote the result of this replacement by 52VF(0). 
(Since both M.2 and M.2 are a factor e smaller than the effective photon mass, 
52W{(t)) and 52W{(t)) actually differ only in subdominant terms.) About the bounce 
solution 52W{(I)) can also be expanded in a derivative expansion, with the leading 
term being the contribution of graph 7b to Z{(j)). 

The third term, lndetG22(0), is the contribution from one- loop graphs with 
only internal 02 propagators. As with 52W , these propagators have the mass M2- 
Summing over graphs with photon loops attached to the 02-line again replaces this 
mass by M2 and replaces G22(0) by 

(522(0) = -n + Mi(0) (5.14) 
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Two points should be noted here. First, the effects of the change of mass are 
not subdominant, as they were in the previous case. Second, this term cannot be 
expanded in a derivative expansion about the bounce. 

Doing the path integral about the bounce and the pure false vacuum and then 
proceeding as in the previous example leads to 



r = e-^° e~^i e~^'' 



detV^'(0b 



detH/^'(0fv) 



-1/2 



/dctG^sM 



-1/2 o2 

^2i^ + ---) (5-15) 



where 



Bo = WoM - W^o(0fv) 



(5.16) 



Bi = 6iWM + d2WM + tr[W^^'(0b)]-^5iW^"(0b) 

(5.17) 

+ trK'(0b)]-'52W^"(0b) - (0b ^ 0fv) 

and B2, containing higher order terms, is of order unity. The leading term in 
the exponent, Bq, contains the part of the effective potential and the tree level 
kinetic term. The four terms in Bi reproduce, in the order in which they appear, 
graphs 7a and 7b of the O(e^) part of and graphs 5a and 5b of the O(e^) part of 
V^eff- Although some of the O(e^) terms in V^g are contained in B2, the potentially 
complex ones, arising from scalar loops, are contained entirely within the explicit 
determinant factors. Because these determinants cannot be expanded in derivative 
expansions, the imaginary part of the effective potential does not explicitly enter 
the calculation. Just as before, it is necessary to look for negative eigenvalues in 
the spectra of the operators in these determinants. By viewing the bounce as the 
optimum tunneling path in the space of field configurations involving both 0i and 
02, it is clear that there is still just a single negative mode, involving only 0i. 
However, in addition to the four translational modes there is a fifth zero-frequency 
mode, corresponding to a phase rotation of the bounce. This must be treated 
by introducing a collective coordinate, replacing detG'22(0b) by det'G22(0b)) 
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multiplying by a Jacobean factor 



J, 



phase 



27r / d^x(^\{x) 



(5.18) 



There remains a minor point, which was noted above. Because the photon 
mass M(0) is proportional to 0, there can be infrared problems if <C cr, as is the 
case far from the center of the bounce, where (\>\,{x) exponentially approaches the 
false vacuum 0fy = 0. These can be seen most clearly by recalling the estimate 
(4.20) for the effective action contribution containing j derivatives and k factors 
of 0. With M(0) = 60, the previous estimate is replaced by 



4+i 



i-4 



(5.19) 



For small 0, the derivative expansion cannot be carried out beyond the four- 
derivative terms. However, if we extract from W the potential terms and the 
two-derivative terms, the remainder is 0(5'^) and gives an 0(1) contribution to the 
exponent S, which can eventually be absorbed in a redefinition of the prefactor. 

To summarize, the bubble nucleation rate is again given by an expression of 
the form of Eq. (4.27), with 



-I 



Kff(0b) + 2(^M0b)' 



(0b ^ 0fv) !> = 0(e-4) 



(5.20) 



and 



Jd^x{ 



K'ff(0b) + ^^;'(0b)(5M0b)' 



(0b ^ 0fv) ^ = 0(e-2) (5.21) 



and the pre-exponential factor A now being of order a^/e^. 
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6. Concluding Remarks 



In this paper I have shown how the decay rate of a metastable vacuum can 
be calculated in a theory whose vacuum structure is determined by radiative cor- 
rections. As in the standard case, the result may be written as a dimensionful 
prefactor times the exponential of an action involving a bounce solution. To lead- 
ing approximation this exponent is just the tree-level action supplemented by the 
dominant one-loop contribution to the effective potential; in scalar electrodynamics 
it is 0(l/e^). The first correction to the exponent arises from the next- to-leading 
contributions to the effective potential and the leading correction to the tree-level 
kinetic part of the effective action. Although smaller than the leading terms, these 
give an addition to the exponent which is larger than order unity (e.g., it is 0(l/e^) 
in scalar electrodynamics) and is thus more important than the prefactor. It does 
not appear that this correction need have any particular sign, but rather that it 
might increase the nucleation rate in some theories and reduce the rate in oth- 
ers. Further, the separation between the leading and next- to-leading terms in the 
bounce action is dependent on the precise specification of the renormalization con- 
ditions. In fact, the existence of a correction term of this magnitude follows simply 
from the requirement that the physical nucleation rate be independent of the renor- 
malization scheme and of the particular definition of the coupling constants of the 
theory. 

All further corrections may be absorbed into the prefactor. Although some of 
these (those corresponding to graphs with internal lines of the heavy "x-particles" ) 
can be identified with particular terms in the effective potential and the other 
functions entering the effective action, this is not true of all the higher corrections. 
Specifically, the 0-loop graphs which give rise to complex terms in the effective 
potential cannot, when calculated in the background of the bounce, be expanded in 
a derivative expansion. Consequently, the imaginary part of the effective potential 
does not explicitly enter the bubble nucleation calculation, and the problems of 
interpretation which it would entail are avoided. In particular, the expression 
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obtained for the nucleation rate can be shown to be real. 

In the examples considered in this paper, it was obvious that the standard 
formalism for vacuum decay had to be modified, since the tree-level potential had 
only a single vacuum and so could not possibly lead to a bounce solution. However, 
one could easily construct examples in which the tree-level potential had several 
inequivalent minima, but where the one-loop corrections to the effective potential 
were comparable to the tree-level terms and changed the vacuum structure. In 
such cases a bounce solution to the tree-level action would exist and could be used 
to calculate the nucleation rate, but the result would not be correct.*^ Working 
with the standard formalism, would one see any indication that the calculation 
was unreliable? The answer is yes. To see this, consider a model similar to that of 
Sec. 4, but with an additional 0^ interaction of magnitude such that F (0, x = 0) 
has two minima. Further, assume that the self-interactions can be written in the 
form 

y(0,X = O)=Aa4t/(0/a) (6.1) 

with U involving no small couplings and the relation between A and the couplings 
involving x being the same as in Sec. 4. Applying the standard formalism in 
this case, one will obtain a bounce solution with nontrivial but with x = 
everywhere and will be led to an expression for the nucleation rate which differs 
from Eq. (2.13) only by a multiplicative factor 



(6.2) 



det[-n + M2(0f^)] 

arising from the functional integration over x- 

The magnitudes of the various terms entering the expression for F can then 
be estimated by the methods of Sec. 2. In particular, since the spatial extent of 



6 A similar situation could arise in the study of solitons. The calculation of the quantum 

corrections to the soliton energy has many similarities to the decay rate calculation consid- 
ered in this paper. Here too, the standard calculation should break down if the one-loop 
corrections to the effective potential are large enough. 
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the bounce solution is of order l/{-\f\a), it is useful to define the dimensionless 
variable s = \/~\ax, in terms of which the bounce has spatial extent of order unity. 



The bounce action is then seen to be of order 1/A, while the determinants arising 
from the functional integration over combine with the Jacobean to give a factor 
of order unity. However, the remaining factor, K^, now becomes 



Because the couplings have been assumed to be large compared to the ^ self- 
interactions, M^(0)/A(7^ is large (of order l/g^). As a result, is not simply 
a factor of order unity which can be included in the pref actor. Instead, it gives 
potentially large corrections whose evaluation'^"' is nontrivial; in other words, the 
standard algorithm has broken down. 

One final note. This work has been concerned solely with the problem of 
bubble nucleation by quantum mechanical tunneling at zero temperature. Issues 
similar to those encountered here also arise in connection with the problem of finite 
temperature bubble nucleation, even in theories where radiative corrections have 
little effect on the zero-temperature vacuum structure. These will be considered 
elsewhere. 
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det[-D, + M2(0b(5))/Aa2] 
det[-n, + M2(0f^)/A(72] 



1-1/2 




(6.3) 
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FIGURE CAPTIONS 



1) Examples of graphs, which although multiloop when drawn in terms of el- 
ementary propagators, are included in the one 0-loop contribution to the 
effective potential in the model of Sec. 4. Solid and wiggly lines represent 
and X propagators, respectively. 

2) The two-loop graphs which contribute to the effective potential in the model 
of Sec. 4 at order fg^ and g^. 

3) The one-loop counterterm graphs which contribute to the effective potential 
in the model of Sec. 4 at order fg^ and g^; the heavy dots denote counterterm 
insertions. 

4) A one-loop graph from which the leading contributions to Z((f)) can be com- 
puted. 

5) The two-loop graphs which give O(e^) contributions to the effective potential 
of scalar electrodynamics. Solid, dashed, and wiggly lines refer to 0i, 02, and 
photon propagators, respectively. 

6) The one-loop graphs which give O(e^) contributions to Zfj^^cf)) in scalar elec- 
trodynamics. 

7) Multiloop graphs whose contributions must be added to 52VF(0). 
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